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I. INTRODUCTION 

The discovery of charmoniuni-like states by the B- 
factories named XYZ mesons, which could not simply 
be described by the naive quark model, have increased 
the interest for new and more accurate charmonium po- 
tentials trying to accommodate these states. 

The static charmonium potential can be parametrized 
phenomenologically as a coulomb plus linear part to- 
gether with spin-spin, spin-tensor and spin-orbit terms 
as leading spin dependent corrections jlj, From a 
more fundamental point of view, this potential can be 
derived from first principles of QCD in the lattice with 
relativistic corrections classified in powers of the inverse 
quark mass 0]. Recently these potentials have been im- 
proved either by new approaches in lattice QCD ^] or by 
matching the long range part calculated by lattice sim- 
ulations of the full QCD with results from perturbative 
QCD at short range [5]. Despite these improvements the 
short range part of all these potentials still suffers from 
singularities which can be only handled by using ad hoc 
regulators. This triggers an unpleasant short distance 
sensitivity. 

In Ref. Q we have developed an approach which re- 
duces the effects of the ad hoc regulators, treating ex- 
actly the singular contributions of the potential following 
renormalization ideas. Instead of fitting the regulators 
to reproduce the ground state mass we use this mass as 
input parameter of our calculation. The rest of the spec- 
trum is predicted from the orthogonality with the ground 
state. The short range uncertainties are encoded in this 
input and the rest of the spectrum only depends on the 
well established pieces of the potential. 
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We have shown in Ref. Q that this procedure accu- 
rately reproduces the same numerical results as the stan- 
dard constituent quark model with extra regulators for 
the particular case of the J^'-' = 1 states. In this pa- 
per we extend the calculation to the full charmonium and 
bottomonium spectrum showing that for all purposes the 
regulators only account for the ground states of cc and 
hb systems with J^'~^ quantum numbers. Therefore the 
approach constitutes a useful tool to complement more 
fundamental potentials. For completeness, we will also 
compare the results with those of the original model (in- 
cluding regulators), where generally an almost perfect 
agreement is found. Once these short range uncertainties 
have been resolved we try to analyze which of the cur- 
rently existing XYZ states could be identified as purely 
quarkonium states. 

The plan of the paper is as follows. Section |lT] ex- 
plains briefly the features of our constituent quark model 
and describes the renormalization approach extended to 
all possible channels of quantum numbers. Sections IIIII 
and IIVI are devoted to comment the spectrum of char- 
monium and bottomonium sectors within the renormal- 
ization approach. Section Hill includes a comparison with 
the original model. We finish summarizing the work and 
giving some conclusions in Section [V] 



II. RENORMALIZED QUARKONIUM MODEL 

The renormalization procedure we use exploits the lo- 
cal character of the potentials and uses a radial regulator 
in terms of boundary condition at a given cut-off radius 
Tc which is made smaller than any other length scale 
of the problem (typically Tj, ^ 0.01 is enough). Within 
the present context the procedure has been explained in 
detail in Q for which we refer for further aspects and 
motivation. Here, we will only present its generalization 
for any channel different than the J^'~^ — 1 . 

The inter-quark potential we use is based on the one 
developed by Vijande et al. which is able to describe 
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meson phenomenology from the hght to the heavy quark 
sector. It also successfully describes hadron phenomenol- 
ogy and hadronic reactions [8l-[lTj|. The quark model 
includes a screened linear confinement potential based 
on unquenched lattice calculations together with spin- 
dependent terms determined by perturbative one-gluon- 
exchange as a Fermi-Breit type interaction. Although a 
complete description of the potential and the parameter 
values are given in Ref. 12\ , it is instructive to write down 
the explicit expressions of the (unregularized) One Gluon 
Exchange (OGE) contribution to display its singularities 
explicitly. 

The OGE potentials are 
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As we see, the tensor and spin-orbit potentials present 
singularities at short distances. Note that we 
have discarded the Dirac delta function (which has 
the same dimensions as which is traditionally in- 

cluded. While this may seem weird, these are distribu- 
tions around the origin which are not seen by the com- 
pact support test functions implied by the boundary con- 
dition regularization below the radial cut-off radius @. 
This also applies to any derivatives of the Dirac delta 
function. The result was suggested jlS] and explicitly 
checked by using a momentum space regularization with 
so-called counterterms [T^ . 

Once the constituent quark model we use has been pre- 
sented, we will apply the renormalization with boundary 
conditions to eliminate the regulators of the model and 
treat exactly the singular contribution of the potential. 



This scheme has been explained in detail in [6| and we 
will provide here only the most relevant additional as- 
pects to deal with all J^^ quarkonium states. 

In Ref. [6*] it is shown that for unique and finite normal- 
izable solutions the number of free independent param- 
eters of the regularized theory depends on the behavior 
of the solution at the origin and can be established by 
a simple analysis of the potentials. The different cases 
that we can find are 

• Uncoupled channels with a singular attractive po- 
tential. In this case, besides the normalization con- 
dition which eliminates one constant, one more pa- 
rameter is needed to renormalize the solution. If 
the potential is repulsive at short range all the con- 
stants which determine the wave function are de- 
fined and the bound state is predicted as it usu- 
ally happens for the standard non singular quan- 
tum mechanics problem. 

• In the case of coupled channels the number of free 
parameters depends of the values of the potential 
near the origin, and more specifically on the corre- 
sponding eigenvalues of the coupled channel poten- 
tial matrix. For two attractive potential eigenval- 
ues we need three observables to fix the wave func- 
tion. If there is one attractive eigenvalue only one 
parameter is needed. Finally, in the case of two re- 
pulsive eigenvalues the wave function is completely 
determined without any additional input parame- 
ter. 

Focusing on quarkonium systems, we will work within 
the nonrelativistic framework, so the dynamics of the sys- 
tem is given by the Schrodinger equation. For a tensor 
and spin-orbit interaction we find the following situations 
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Triplet uncoupled channel {s = 1, I = J) 
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Triplet coupled channel (s = l,^ = J±l) 
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where the energy is defined with respect to the Q — Q threshold, £"„ 



Mn - 2m. 
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As mentioned, in order to determine the number of 
independent constants we have to study the potential at 
short distances for different channels. In the r — )■ limit 
the dominant contributions are the tensor and the spin- 
orbit terms of the OGE potential. We need to know their 
character for the different cases: 

• Singlet channel (s = 0, ^ = J). We have (5i2) = 
and {L ■ S) = 0, so the potential is regular and 
there is no observable free parameter. For the ^Sq 
channel we have to take into account that in this 
case the potential has an attractive (5-function and 
so is singular, therefore we have to set an observable 
for a regularized solution. 

• Triplet uncoupled channel (s = 1, Z = J). We have 
('S'12) = +2 and {L ■ S) = —1, so the potential is 
singular attractive and an observable must be fixed 
for a regularized solution. 

• Triplet coupled channel (s = l, Z = J±1). If we 
denote / = J — 1 and Z' = J + 1 , we will have 



- eij\L-s\^ij 



{%\L-S\%)^-iJ + 2), 

and diagonalizing the potential matrix the eigen- 
values are 



E{J) = -10 ± 6 VT + 7+J2 



(7) 



with J > 1, so we always have one negative eigen- 
value and it requires one parameter. 

In order to describe a bound state we seek normalizable 
solutions 



(8) 



where w{r) = for the uncoupled channels. This imposes 
conditions on the wave functions both at infinity as well 
as at the origin. 

Moreover, the set of equations dH [SJ [H]) must be ac- 
companied by asymptotic conditions at infinity. Once 
we have discard the irregular function at long distances, 
the wave functions at infinity have the following behavior 



u{r) — Aj^i e 
w{r) Aj+i e~ 



(9) 



where Aj^i is the normalization factor and the asymp- 
totic J + 1/ J — 1 ratio parameter is defined by 7y = 
Aj+i/Aj-i. Ideally, one would integrate the Schrodinger 
equation taking its solutions at infinity, Eq. which 



depend on the binding energy and r]. The singular struc- 
ture of the problem at short distances requires a specific 
analysis of the coupled equations as it has been done 
extensively elsewhere [15| and we adapt here for our par- 
ticular situation. The result amounts to integrate from 
infinity for the physical value of A/q and ry (or Mq in the 
case of singlet channels). Generally, the solutions diverge 
strongly at the origin, so that the normalization of the 
state is precluded. However, there is a particular value 
of 77 which guarantees that the wave function becomes 
normalizable ^. The rest of the spectrum is then built 
by imposing orthogonality of states in coupled channels 
(see i§] for details in the J ftp case). 



III. PHENOMENOLOGY OF CHARMONIUM 
STATES 

The charmonium spectrum consists of eight narrow 
states below the open-charm threshold (3.73 GeV) and 
several tens of states above the threshold, some of them 
wide, because they decay into charmed mesons, some of 
them still narrow, because their decay to open-charm is 
forbidden by some conservation rule. Below the thresh- 
old all states are well established. Above threshold, how- 
ever there are new charmonium-like states that are very 
difficult to accommodate theoretically. 

In Tables 4 and 9 of Ref. JJij the updated new con- 
ventional and unconventional states in the cc, be and bb 
sectors are given. If we focus on the cc region, there 
are three new states which have been recognized as qq 
pairs in the PDG [l3] during the last years. They are the 
he which is the ^Pi state of charmonium, singlet part- 
ner of the long-known Xcj triplet ^Pj states, the r]c{2S) 
which is the first radial excitation of the pseudoscalar 
ground state ric{lS) and the Z(3930) whose assignment 
as the 2^P2 state, Xc2(2P), seems widely accepted. The 
rest of the resonances of this sector, namely the X{3872), 
X(3915), r(3940), A:(3940) stiU lack for a clear assign- 
ment. 

The results of the study of the full charmonium spec- 
trum up to total spin J = 2 and for the first radial exci- 
tations within the renormalization with boundary condi- 
tions scheme are shown in Table HI In the following we will 
discuss our predictions for the different channels except 
the J^'~^ = 1 which has been extensively discussed in 
Ref. already. 

In columns four and five we compare the calculated 
masses within the renormalization approach with the ex- 
perimental data. As we see, once the experimental value 



^ Thus, if one imposes the regularity condition at the origin one 
will determine 77 and therefore the wave function of the bound 
state. In practice, however, the converging solution is rather 
elusive since integrated-in solutions quickly run into diverging 
solution due to the round-off errors and dominate over the con- 
verging solution. 
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of the ground state mass is taken as input in the calcu- 
lation, the rest of the spectrum is reproduced fairly well. 

In order to size the relevance of form factors in the 
excited spectrum we also show in columns six and seven 
the comparison with the original model including form 
factors [7\ and taking the predicted mass as a parameter 
of the renormalization approach. The observed tiny devi- 
ations are below a few MeV corresponding to a marginal 
influence of the form factors on the excited states. These 
results actually point to the idea that the gluonic regula- 
tors are fitted just to provide the ground state energies. 
Once we have eliminated possibles bias due to the use 
of regulators, we are in a favorable position to discuss 
possible qq assignments for the old and new states. 



A. rjc and ric2 states 

An ric{lS) candidate has been observed thirty years 
ago by CBAL 21] and MARK II ^ CoUaborations with 
a mass measurement very close to the updated world av- 
erage 2980.3 ± 1.2. In the renormalization scheme this 
state is taken as a parameter and in essence fixes the 
spin-spin contact term interaction. 

The search for a reproducible rye (25) signal has a 
long history. Recently, Belle [l^ found a signal in 
B -J> Kr]c{2S) in the exclusive r]c{2S) K°gK-T:+ de- 
cay mode (a favorite all-charged final state for 7/c(l<5')), 
at 3654 ± 6 ± 8 MeV. Since then measurements of r]c(2S) 
in that mass region have been reported by BaBar [24| . 
CLEO [ii], and Belle 26] in 77-fusion to KKtt final 
states and by BaBar [2^ and Belle ^28] in double char- 
monium production. 

Our predicted mass for the r]c{2S) is 3634 MeV in very 
good agreement with the updated world average reported 
in Ref. [l^. Nothing is known of the next excitation, 
rjci'iS). Our prediction is around 4.05 GeV. 

The potential at short distances is regular for the ric2 
states. Therefore the mass of the ground state for the r]c2 
meson is not a parameter and we predict 3812, 4166 and 
4437 MeV for the ground state and the first two radial 
excitations. 

We can clearly identify a second state ric2{2S) with 
mass M = 4166 MeV and width T = 122.9 MeV with the 
resonance recently reported by Belle at M = 4156^20 ^ 
15 MeV with a width T = 139^11^ ± 21 MeV [29^ in the 
e+e~ — ^ D*D*J/Tp. The decay of r]c2 to DD is forbid- 
den being the X(4160) D*D* the most favored decay 
channel as shown by the data. 



B. he and Xd states 

Two experiments reported the observation of the 
hc{lP) in 2005. CLEO [H HH obtained a 6a statis- 
tical significance in the isospin-forbidden decay chain 
e+e" tP{2S) ^ 7r°/ic, iVci^S). E835 ^ found 

a 3(7 evidence in pp — he, he — jrid^S), r]c{lS) — 77. 
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TABLE I: Masses, in MeV, of charmonium states. nL labels 
the radial and angular momentum quantum numbers. When 
partial waves are coupled they refer to the dominant compo- 
nent. We compare with the well established states in Ref. [13] 
and assign possible XYZ mesons. 



The precision measurement of its mass was reported 
by CLEO in 2008 [H, 3525.28 ±0.19 ±0.12) MeV. Later 
BES III "34^ has confirmed this with a mass of 3525.40 ± 
0.13 ± 0.18 MeV. It was important to measure the mass 
of this state because the spin-averaged ccntroid of the 
triplet states 
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TABLE II: The theoretical masses, in MeV, of the first two 
radial excitations of he compared with the spin-averaged cen- 
troid, in MeV, of the triplet states. We compare with the 
experimental data [l7j . 



is expected to be near the hc{lP) mass. The lattice 
data show a vanishing long-range component of the spin- 
spin potential. Thus, the potential appears to be entirely 
dominated by its short-range, delta-like, part, suggesting 
that the ^Pi should be close to the center-of-gravity of 
the "^Pj system. So it makes the hyperfine mass splitting, 
Amhf[hc{lP)] = {m(l^Pj)) - m[hc{lP)] an important 
measurement of the spin-spin interaction. 

The centroid of the l^Pj states (xo.i.s) is known to 
be [H] 3525.30 ±0.04 MeV and then the hyperfine split- 
tings are ±0.02 ± 0.23 MeV, from CLEO and -0.10 ± 
0.22 MeV from BES IIL 

Table U shows the masses for three radial excitations 
of the singlet ^Pi and the triplet •^Pj mesons. In Ta- 
ble |TT] we show the comparison between the centroid of 
■^Pj states and the corresponding he mass for the ground 
state and the first radial excitation, showing that our 
spin-spin interaction is negligible for these channels and 
it is in perfect agreement with the lattice expectations 
and the experimental measurements for the ground state. 

The mean 2P multiplet mass is predicted to be near 
3.95 GeV. Although no 2P cc states have been clearly 
seen experimentally, there are a lot of states reported 
from the different Collaborations which claim enhance- 
ments in that energy region. Concerning the rest of the 
new states, none of them appears in our calculation as a 
cc state except the ^(3930). This fact agrees with the re- 
sults of a recent calculation which shows that these states 
can be interpreted as molecular resonances [ssj . 

The Z(3930) was reported by Belle in 77 -> DD with 
a mass and width M = 3929 ± 5 ± 2 MeV and T = 
29.9±10±2 MeV 'll|. The two photon width is measured 
to be T^^B{Z {3930) DD) = 0.18 ± 0.05 ± 0.03 keV. 
Moreover the DD angular distribution is consistent with 
J = 2. The Xc2(2P) state is a good candidate for the 
Z(3930). We obtained a mass of 3974 MeV and the total 
width r — 49.1 MeV is comparable with the experimental 
data. Finally the experimental two photon width com- 
pares nicely with our result T-f^B{Z{3930) — > DD) = 
0.15 keV. 



IV. BOTTOMONIUM 

In this Section we extend the renormalization with 
boundary conditions scheme to the study of the bottomo- 
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TABLE III: Masses in MeV of bottomonium states. The no- 
tation for the quantum numbers is the same as in Table [iTl 
We compare with the well established states in Ref. [TtI ]. 



nium states. Although the S-factories are not usually 
considered as ideal facilities for the study of the bottomo- 
nium spectrum since their energy is tuned to the peak of 
the T(4S') resonance, which decays in almost 100% of 
cases to a BB pair, BaBar and Belle have reported data 
samples at various energies in the bottomonium region 
that made possible discoveries like the r^f, [33|, ht,{lP) 
and hi,{2P) Nevertheless the number of states is 
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TABLE IV: The theoretical masses, in MeV, of the first two 
radial excitations of hb compared with the spin-averaged cen- 
troid, in MeV, of the triplet states. We compare with the 
experimental data [l7j . 



in reasonable agreement with the experimental data. In 
Table ITVl we show the comparison between the centroid of 
Xbj states and the corresponding hb mass for the ground 
state and the first radial excitation, showing that our 
spin-spin interaction is negligible. 

Xf,j(lP) and Yftj(2P) with J = 2, 1 were discov- 
ered earlier in 1982 ^ and 1983 [IH, EH, respec- 
tively. Their masses have not changed much since then 
and our theoretical prediction through both schemes are 
very close to the experimental values. 



lower than in the charmonium spectrum. 

As in the previous Section, column four of Table ITTTl 
shows the calculated masses within the renormalization 
approach up to total spin J = 2 and for the first radial 
excitations. As in the case of the charmonium spectrum, 
when applicable, we indicate the fact that the ground 
state mass is used as an input parameter. 

In the J^^ = 1^^ channel there are less data for the 
bottomonium spectrum than for the charmonium. The 
S and D wave states are almost degenerate for the high- 
est excited states. In the charmonium these doublets 
has been measured and we can assign them to S and 
D waves. In the bottomonium the splittings are smaller 
and experimentally they have not been resolved. We fol- 
low the assignment of the Particle Data Group up to 45' 
states and assumed that the others are S wave. This is 
the reason why we assigned the experimentally measured 
state at M = 11019 MeV to our state 65 and not to our 
5D. 

As a general trend the experimental data are well re- 
produced by the calculation and a number of excited 
states are given which can be tested at LHCb and other 
B-factories. The rii,(2S) has been very recently mea- 
sured by the Belle Collaboration [sBl with a mass M — 
9999.0 ± 3.5tl^ which is in reasonable agreement with 
our M = 9957 MeV resuh. 

The hyperfine mass-splitting of singlet-triplet states, 
Amhf[rib{lS)] = m{T{lS)) — m{rib{lS)), probes the spin- 
dependence of bound-state energy levels, and, once mea- 
sured, imposes constraints on theoretical descriptions. It 
is given experimentally by 

Am,,/ [ryb (15)] = 69.6 ±2.9 MeV (11) 

In the renormalization scheme this splitting is fixed 
and we can only predict the splittings for radial excita- 
tions which has not yet been measured. 

In the case of the centroid of the Xbj{nP) states 
with n = 1, 2 is known to be [13) 9899.87 ± 0.27 MeV 
and 10260.24 ± 0.36 MeV, respectively. The hyrcrfine 
splittings measured by the Belle Collaboration [38] are 
Am,,/ [/lb (IP)] = +1.6 ± 1.5 MeV and Am,,/[/ib(2P)] = 
+0.5l|^'2MeV which are compatible with zero. 

Table IIIII shows the masses for three radial excitations 
of the singlet hb and the triplet Xbj mesons. They are 



V. CONCLUSIONS 

Based on an earlier work, we have reanalyzed the cal- 
culation of the charmonium spectrum in a constituent 
quark model using a renormalization scheme with bound- 
ary conditions. This approach avoids explicitly the in- 
troduction of phenomenological form factors taking as 
a parameter the mass of the ground state. Thus, the 
only relevant physical information on the form factors 
is to tune the value of the ground state energy. Once 
this fact has been established we have applied the renor- 
malization framework to provide some basis to "bare" qq 
assignments of mesonic states. 

We obtain a spectrum in reasonable agreement with 
the experimentally well established data. For instance, 
we obtain Am,, f [hc{nP)] compatible with zero due to the 
fact that the hyperfine contact term has been included 
in the renormalization conditions. We also assign certain 
XYZ mesons according to our model. 

For the phenomenologically successful model of Ref. 
where ad hoc form factors are introduced as regulators, 
we find an almost perfect agreement with the renormal- 
ization approach. This result provides confidence on the 
way the original model took into account the unknown 
short distance dynamics. In addition, we have extended 
this study to the bottomonium sector obtaining similar 
conclusions as in the charmonium sector. 
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